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Abstract 

We compare the behaviour of a new 2D aqueous colloidal model system with a simple numerical 
treatment. To first order the attractive interaction between the colloids induced by an in-plane 
rotating AC electric field is dipolar, while the charge stablisation leads to a shorter ranged, Yukawa- 
like repulsion. In the crystal-like 'rafts' formed at sufficient field strengths, we find quantitative 
agreement between experiment and Monte Carlo simulation, except in the case of strongly interact- 
ing systems, where the well depth of the effective potential exceeds 250 times the thermal energy. 
The 'lattice constant' of the crystal-like raft is located approximately at the minimum of the effec- 
tive potential, resulting from the sum of the Yukawa and dipolar interactions. The experimental 
system has display applications, owing to the possibility of tuning the lattice spacing with the ex- 
ternal electric field. Limitations in the applied field strength and relative range of the electrostatic 
interactions of the particles results in a reduction of tunable lattice spacing for small and large 
particles respectively. The optimal particle size for maximising the lattice spacing tunability was 
found to be around 1000 nm. 

PACS numbers: 82.70.Dd 
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I. INTRODUCTION 



Colloidal suspensions present the possibility to develop novel materials via self-assembly. 
Of particular interest are colloidal crystals, whose optical properties can generate iridescent 
colours, and provide a means by which photonic crystals may be produced(author?)-, 
while further applications range from lasers (author?) 2 to display devices (author?)-, 
with recent advances demonstrating tunable colours through control of lattice spacing with 
an external field (author?) 4 . Further to the practical importance of colloidal crystals, 
their well-defined thermodynamic temperature allows colloidal dispersions to be viewed as 
mesoscopic 'model atoms' (author?) 5 . 

Recently, the ability to tune the colloid-colloid interactions has led to the observation of 
a wide variety of structures^ 1 ^^^. Of particular interest here, to first order AC electric 
fields can induce dipolar interactions between the colloidal particles, leading to anisotropic 
interparticle potentials and exotic crystal structures, some of which are not observed in 
atomic and molecular systems (author?)-, while external control of the colloid-colloid 
interactions allows direct observation of phase transitions^. Furthermore, direct micro- 
scopic observation at the single-particle level allows an unprecedented level of detail to be 
accessed^, opening the possibility of tackling long-standing problems in condensed matter, 
such as freezing^. 

The introduction of a rotating AC field opens up even more possibilities. In this case, the 
dipolar interactions lead to an attraction in the plane of rotation and to repulsions above 
and below. Studies with a rotating magnetic field on granular matter indeed produced 
disc like patterns consistent with expectations^ 5 - 1 ^. Unlike granular matter, since colloidal 
dispersions exhibit Brownian motion, thermodynamic equilibrium structures (ie crystals), 
may be obtained (author?)—'—. In previous work Snoswell et.al. (author?)— showed 
that lattice spacing within quasi-2D colloidal crystals could be controlled in-situ, by means 
of coplanar rotating electric field. The interparticle dipolar interactions in the plane of the 
electric field may be treated to first order as a circularly symmetric attraction, due to the 
time averaging effect of a rapidly rotating field (1000Hz) on relatively large particles on the 
micron lengthscale, where the diffusive timescale is of the order of seconds (author?) 18 . 

In considering the interactions between the particles, the asymmetry between the col- 
loids (10nm-l/xm) and smaller molecular and ionic species must be addressed. A number 
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of coarse-graining schemes have been developed where the smaller components are formally 
integrated out (author?) 19 . This generates a one-component picture, where only the ef- 
fective colloid-colloid interactions are considered, and the complexity of the description is 
vastly reduced. The equilibrium behaviour of the colloids in the original multi-component 
system may then be faithfully reproduced by appeal to liquid state theory (author?)— and 
computer simulation (author?) 21 . Central to the success of this one-component approach 
is the use of a suitable effective colloid-colloid interaction u(r). 

In this study, we use a simple numerical treatement in which we can predict the lattice 
spacing in the quasi 2D crystal from the electric field strength. We consider a model system 
of charged colloids, in a rotating electric field (author?) 11 . By exploiting the knowledge 
both of the electrostatic repulsions and dipolar attractions, we present a direct, quantitative 
comparison of a tunable interaction and a material property of the crystalline 'rafts' formed. 
We combine experimental measurements of the crystal lattice constant d as a function of 
field strength E with Monte-Carlo simulations according to a screened Coulomb repulsion 
plus dipolar attraction where the only fitting parameter is the Debye screening length. 

In the simulations, we use pairwise interactions, in other words we assume that at the 
higher densities at which the crystalline rafts are formed, the system is still accurately 
described by interactions calculated for two particles in isolation. We note that deviations 
from this assumption of pairwise additivity have been measured both in the case of strongly 
charged colloids (author?)— and in the case of repulsive dipolar interactions^. We further 
compare simulation results with the minimum of the effective potential, which we take as 
a measure of the lattice constant of the crystalline rafts, which we also determine from 
experimental data. 

This paper is organised into six sections. In section [Til we present expressions for the 
effective interactions between the colloids, summing the attractions and repulsions to provide 
an effective one-component description of the system. Section ITlTI describes our experimental 
metholodogy. Section IIVI outlines the Monte-Carlo simulation technique employed. The 
comparison of simulation and experimental results is presented in section |V] and in section 
IVII we extrapolate our findings to maximise the tunability of the crystal lattice constant, 
which may be useful for applications. We conclude our findings in section IVHI 
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II. THEORY AND MODEL 



In the following we will consider a system consisting of two particles in a surrounding 
medium. We shall assume that these particles are charged, leading to a repulsive interac- 
tion, and that the rotating AC electric field induces a dipole moment in the two particles 
and thus induces an attractive interaction. To describe this system, we start from the Der- 
jaguin, Landau, Verwey and Overbeek (DLVO) approach 24 , which consists of attractive van 
der Waals interactions at short range, and long-ranged repulsive electrostatic interactions. 
The van der Waals interactions are very short-ranged, and are neglected, as electrostatic re- 
pulsions inhibit the close approach at which van der Waals interactions become important. 
We shall therefore assume that the only relevant attractions result from the long-ranged 
dipolar interactions induced by the rotating electric field. 

In the linear Poisson-Boltzmann regime, the electrostatic repulsions may be expressed as 
a hard core Yukawa, or screened Coulomb interaction 2 ^, 



where (3 = l/ksT where ks is Boltzmann's constant, T is temperature and a is the colloid 
diameter. The potential at contact, (3e yu k is given by 



where Z is the colloid charge, Ib is the Bjerrum length and k = y/ATrlsPi is the inverse Debye 
screening length where pi is the total number density of monovalent ions. 

Now the regime of linear Poisson-Boltzmann theory in which equation (OQ) holds corre- 
sponds to relatively weak charging. Although this is not the case here, a potential of the 
Yukawa form is recovered at larger separations, if a smaller, renormalized charge is con- 
sidered (author?) 25 . We tabulate measurements of the (^-potential of dilute suspensions 
in table [B These values suggest that we expect a renormalised charge, in the conditions 
under which the colloids form the crystalline rafts, rather high colloid concentration, hence 
high counter ion concentration 1 ^ 5 -. Therefore, noting that the Debye length is much smaller 
than the colloid radius, we follow Bocquet et. alM- and take the following expression for the 
renormalised charge Z e ff : 




for r < a 



for r > a 



(1) 




Z 2 l B 
(l + na/2) 2 a 



(2) 
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4a (1 + Ka/2f 
Zeff ~ T B 1 + na (3) 
which we substitute for Z in equation (j2j). This expression gives good agreement with 

measurements of the effective colloid charge for particles with a comparable C-potential^I. 
We recall that many-body effects can lead to a density- dependence in the effective colloid- 
colloid interactions (author?)—'—. However, we shall neglect these effects in the present 
work. 

The attractive potential between the colloids resulting from the rotating AC electric field 
is treated to first order dipolar interaction: 

(4) 

where p is the induced dipole moment of each particle, e is the dielectric constant and eo is 
the permittivity of free space. The dipole moment may be calculated from the strength of 
the electric field 

p = ~7ree a 3 KE (5) 

where K is the Clausius-Mosotti factor and takes values between 1 and -1/2 depending 
upon the origin of the dipole moment. We note that in this case of an alternating field, 
the root-mean-square of the time dependent field is taken. Substituting r = a leads us to a 
potential at contact 



edipia) = — 3 (6) 



(g X P) 

27ree cr 3 

which has a cubic dependence on the colloid diameter, in the case of all other contributions 
being unchanged. 

In considering the value of the Clausius-Mosotti factor, two regimes are relevant, corre- 
sponding to applied electric fields of high and low frequency, in which the dipolar interactions 
result from dielectric permittivity differences between the particles and the solvent and dif- 
ferences in conductivity respectively^. The crossover frequency u* between these regimes is 
given by 



where s are conductivities and the indices p and m refer to the colloidal particles and the 
medium respectively. Following (author?)— the conductivity of the particles s p is taken to 
be the sum of the bulk conductivity s& ~ and the conductivity on the surface s s . 

s p = s b + s s a (8) 

In principle equation (jSJ) holds for particles larger than 1 /mi. In the case of much smaller 
particles, one should take into account that the contribution of the diffuse layer to the 
conductivity of the particles. Noting that the smallest particles used in the experiments 
are 757 nm in diameter and that the conductivity of the de-ionised water is relatively low, 
neglect the double layer contribution and thus use equation (jSj). Since the largest particles 
studied are 2070 nm in diameter and the frequency is 1 kHz and the bulk ionic strength 
is 0.01 — 0.1 mMol as determined from conductivity measurements, we work in the low 
frequency regime where the contributions from the conductivities dominate (author?)—. 
The Clausius-Mosotti factor therefore takes the following form: 

K - = T^t (9) 

°p z " 5 m 

Due to the surface conductivity, s p is much larger than the conductivity of the medium 
Sfe. The value of K s is therefore close to one. We now combine the contributions from the 
electrostatic repulsions and the dipolar attractions, yielding the expression. 

u t ot(r) = u yuk (r) + u dip (r) (10) 

Utot has the form indicated in Fig. [U We see a minimum in the potential, which one might 
expect to provide a first approximation to the lattice constant in the 2D colloidal crystal. 
In sections |V] and IVII we shall compare this minimum to our simulation results. 



III. EXPERIMENTAL 



Colloidal crystals underwent self-assembly as a result of an applied electric field. We used 
anionic, sulphate stabilised polystyrene latex particles, either synthesised using a standard 
technique, surfactant free emulsion polymerisation (author?)— in the case of a =757 and 
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945 nm, or particles produced by the same method but purchased from Microparticles GmbH 
for a =1390 and 2070 nm. 

Particle electrophoretic mobility was measured using a Brookhaven Zetaplus light scat- 
tering instrument. Particle sizes were determined by scanning electron microscopy, either 
in-house using a Jeol JSM-6330F, in the case of a =757 and 945 nm or by Microparticles 
GmbH for o =1390 and 2070 nm are are listed in table [B Experiments were performed with 
dilute aqueous suspensions (0.5-1.5 wt%). Schematics of the experimental set up are shown 
in Fig. [2] All glass surfaces were chemically washed with 0.1M KOH and washed with copi- 
ous quantities of MilliQ water. Particles were deionised by direct contact with ion exchange 
beads before being made up to the desired electrolyte concentration with KC1. Very low 
salt conditions are required as even moderate (mMol) salt concentrations lead to electrohy- 
drodynamic pattern formation due to ion flow^^i. For further details of the experimental 
set up and procedure the reader is referred to Snoswell et. a/.—. 

We consider our experiment as a 2D system. However, occasionally we noticed some 
overlap of the crystalline rafts. We do not include these results in our analysis. In order to 
treat the system in 2D, one might expect the gravitational length l g = ksT '/mg, where m 
is the bouyant mass of the particle and g is the acceleration due to gravity, to be much less 
than some characteristic length such as the particle diameter. Due to the fairly small density 
mismatch between polystyrene and water, in fact even for the largest particles a = 2070 nm, 
l g ~ 0.87a. However, the dipolar interactions between the particles are attractive in-plane, 
but strongly repulsive in the vertical direction. This promotes the formation of 'rafts' and 
'sheets', and these large assemblies of many particles have very small gravitational lengths. 
Thus we argue that the system behaves in a quasi 2D manner. The lattice parameter d was 
taken as the average of typically ten crystalline 'rafts', measured across the 'raft', of around 
ten lattice spacings. The response to changing the electric field strength was determined to 
be less than 100 ms, we waited 5s after changing the field strength before acquiring data. 

IV. MONTE-CARLO SIMULATION 

We use standard Monte-Carlo (MC) simulations in the NVT ensemble^. The particles 
interact via equation ( fTOi) in two dimensions. To mimic the experiments, we initialise the 
system in a random configuration at a relatively low concentration, corresponding to an area 
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fraction <p = 0.05 — 0.3. We confirmed that different area fractions gave indistinguishable 
results, however smaller area fractions often led to longer equilibration times. The attractive 
interactions cause the particles to approach one another, and form crystallites which then 
coalesce to form a crystalline 'rafts' which contain of order 100 particles, in a qualitatively 
similar way to the experimental system. Each simulation was equilibrated for typically 
30000 MC moves per particle, followed by 3000 production moves per particle. We recall 
that a potential of the form 1/r 3 converges in 2D. The lattice constant for this 2D hexagonal 
crystal was taken as 

d = ! { r9ir)dr (11) 
Jo 9(r)dr 

where g(r) is the pair correlation function, and b the minimum between the first and second 
peaks of g(r). 

We found that when we used N < 200, that typically one crystalline domain was formed, 
around the centre of the simulation box. We therefore do not use periodic boundary condi- 
tions, and consider the lattice spacing of the single crystal formed in the simulation. One 
parameter is the number of particles typically present in each crystalline domain. This is 
estimated to be around N = 100±25 from experimental data (see Fig. [3]). This value is 
governed by the overall concentration of the colloidal suspension. We have considered the 
effects of varying N as shown in Fig. HI The results of a considerable range of N resulted 
in a variation of a around one percent in d for the a = 945 nm system for a field £"=20 
kVm~ l . Each simulation was repeated four times except iV = 576. We found a slight trend 
to tighter binding for larger N, however, as shown in Fig. HI this effect is rather small. We 
henceforth use one simulation per state point unless otherwise stated. The slight scatter 
in the simulation data is perhaps indicative the existence of different 'magic numbers' for 
these crystalline rafts 32 . These may be thought of as 2D hexagonal clusters, whose 'magic 
numbers', ie low energy states, are expected to include 7, 19, 37... Close to a 'magic number' 
the binding may be expected to be relatively tighter. A more detailed exploration of this 
phenomenon lies beyond the scope of this work. Overall, the variation in d as a function of 
N is smaller than the experimental scatter, nonetheless we take both N = 72 and N = 144 
values for d when comparing with experimental data. The slightly larger value for N = 32 in 
Fig. H]is attributed to the small size of the crystalline 'raft', such that the surface particles 
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which are more widely spaced make a greater contribution. The domain size is of order 100 
particles, thus, for larger N, defects and grain boundaries may lead to a smaller contribution 
to the measurement of d. This also applies to experimental data. 

For fitting, all parameters are known, except the Debye length k _1 which we take as a free 
parameter for each particle size and salt concentration, although we note that the effective 
colloid charge is itself a function of the Debye length, equation ([3]). It has previously been 
shown that comparison of simulation and experimental data can yield a reasonable local 
measure of the Debye length^. Furthermore, in the region of the sample in which the 
electric field is applied, the colloid volume fraction is much higher, and due to the colloidal 
counter-ions the ionic strength may increase, leading to a reduction in the Debye length with 
respect to the bulk. However, it is the Debye length in this region which is relevant for the 
effective colloid-colloid interactions. The bulk Debye length that may be determined, for 
example from conductivity measurements, may therefore be taken only as an upper bound. 
System parameters for different particle sizes are tabulated in table [H and the Debye length 
is plotted in Fig. [5](a) and the contact potential of the Yukawa interaction [equation (T5])] in 
Fig. EJ^b) are plotted as a function of particle diameter. 

Clearly, our MC simulations may be expected to provide a reasonable treatment of the 
crystalline rafts in (quasi) equilibrium, rather than to describe the formation process. We 
therefore restrict our analysis to a simple characterisation of the crystal rafts, and primarily 
consider the lattice spacing. We furthermore note both here and in previous work— that 
considerable variation in shapes of crystalline rafts was observed, but that this had little 
impact upon the measurement of the lattice spacing d. Likewise, at the level of this work, 
we neglect possible local variations in d due to the proximity of an interface. In any case 
we note that for display applications, all lattice spacings contribute to the diffraction. A 
few words on equilibration are in order. This applies both to the experimental system, 
and to the simulations. Neither system is strictly in equilibrium, in that case we might 
expect a rather regularly-shaped raft such as a hexagon. However, the insensitivity of the 
lattice parameter either as a function of N (Fig. H|), and the very close agreement between 
statistically independent simulation runs lead us to conclude that our approach is sufficient 
to compare lattice parameters between experiment and simulation. 



9 



V. RESULTS 



The effect of changing the electric field strength is readily demonstrated in Fig. [3j which 
shows optical microscope images of the region of the sample in which the field is applied. 
Two field strengths are shown, 29 kVm -1 (a) and 80 kVm" 1 (b) for a = 945 nm and we see 
a correspondingly tighter lattice in the case of the higher field strength. Similar behaviour 
was observed ir*^. Figure [S^c) illustrates simulated coordinates for e yu k= 5.78 x 10 3 k B T, 
ko = 28.4, Sdi P = 50.1 &bT, which corresponds to E=29.0 Vm _1 for the a = 945 nm system 
(see table [I]). 

A more quantitative comparison between experiment and simulation is shown in Fig. 0, 
which plots the radial distribution function g(r). We see that the structure appears to be 
crystalline, with reasonably long-ranged correlations, although, as is clear from Fig. [3], the 
system is too small to test for truly long-ranged order. We thus note that it is difficult to 
distinguish the true crystal from a hexatic phase in this system. However, our motivation is 
to model the experimental system, with N = 100 ± 25 and thus we argue that identification 
of the hexatic phase lies beyond the scope of this work. 

Having illustrated the general behaviour of the system, we now turn our attention to the 
fitting of the experimental data, with the model described in section II. The experimental 
data was fitted to the theory by taking the Debye-length as the fitting parameter. 

The main results are shown in Figs. [7] and [8] for o = 757 and 945 nm, and a = 1390 and 
2070 nm systems respectively. These concern the lattice parameter d as a function of applied 
field strength E. In general, the simulation is able to capture the behaviour of the system in a 
reasonably quantitative manner. Furthermore, simply plotting the minimum in the potential 
Umin given by equation (fTOl) provides a first approximation to the lattice constant d. Figure 
[7]^a) shows the results for a = 757nm particles. Due to the small particle size, and the 
inverse cubic dependence of the strength of the dipolar attractions on the particle diameter, 
as characterised by e^ p [equation ([6])] the attraction for a given field strength is comparitively 
small, so we find relatively larger lattice parameters for this system. At low field strengths 
(E < 20 kVm" 1 ) we see an increase of the lattice constant for both the experimental data 
and the simulation, compared to the minimum in the potential. Apparently, as the crystal 
approaches melting, fluctuations become more pronounced, leading to an increase in d, which 
is not captured in equation (flOl) . or perhaps a molten surface layer increases the apparent 
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value of d. We leave this intriguing question for further investigations. Meanwhile, at high 
field strengths, we find some deviation between simulation, experiment and equation ffTOl) . 
Apparently, due to the long ranged nature of the interaction, second nearest neighbours 
experience an attraction of sufficient strength that the crystal is compressed by its own 
cohesion, leading to a smaller lattice constant. 

In the case of the a = 757 nm particles, we find that for small field strengths, (E = 10 
kVm -1 ), the system does not form crystallites, rather it remains as a colloidal liquid. We 
determined this both by experimental observation and with simulations. In the latter case, 
we identified crystallisation with a splitting in the second peak of the pair correlation function 
g(r). The kink around E ~ 20 kVm^'m Fig. [7] is likely related an artifact of measuring 
multiple rafts, so defects and grain boundaries contribute to the value of d. At weak field 
strengths the rafts tend to be smaller, thus increasing this effect. 

At larger particle sizes, (Fig. [8]), we again see reasonable agreement between experiment 
and simulation, albeit with some deviation between simulation and equation ( flUj) at higher 
field strengths, as the simulation predicts a smaller lattice constant than that observed in 
the experiments. However, for the o = 1390 nm system in particular, in fact equation (JTUJ) 
provides a more accurate description of the experimental data. Apparently some of our 
assumptions in the simulations, perhaps that of pairwise additivity, begin to break down 
at high field strengths. According to equation fflOl) . the potential at the minimum of the 
attractive well is some 250 ksT for an applied field of 50 kVm -1 for a = 1390 nm. Thus 
we conclude that for moderate interaction strengths, equation ( flOi) provides a reasonable 
description of the system. 

We now consider the dependence of the system upon the colloid diameter a. Firstly, we 
see from Fig. [5]^a) that the absolute value of the (fitted) Debye length k~ 1 does not change 
significantly for all the four particle sizes studied. Thus, the reduced inverse Debye length 
kg is linear in a . The contact potential (3e yu k has an approximately inverse cubic dependence 
on a for kct^> 1, equation (T5]), but also depends upon the (effective) charge Z e ff. Plotting 
(5e yu k as a function of a in Fig. [5](b) we find an approximate power law dependence with an 
exponent of b = 0.91 ± 0.01. 

As noted above, the prefactor of the dipolar attraction, e^p [equation (151)], has a ct~ 3 
dependence. Although the electrostatic interactions are non- negligible, we nonetheless ex- 
pect from equation that upon decreasing a we find a relatively larger lattice constant 
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for a given field strength, and that larger field strengths are required to provide sufficient 
interactions that the crystalline rafts form. This we indeed find, as shown in Figs. [7| and [H 
and also in the next section. 

VI. OPTIMISING THE TUN ABILITY 

Now we have noted that this system has optical display applications, due to the possibility 
of externally tuning the lattice parameter d. These applications may be most usefully realised 
when the possibility to tune the system is maximised, ie when the range of d is maximised. 
We are therefore motivated to consider a range of particle sizes, and calculate the range of 
tunability of d. Now a sets a lower bound to d while the upper bound is set by the melting 
transition which can be determined by Monte-Carlo simulation. 

Hitherto, we have used the Debye length fitting parameter. However, as Fig. 

E^a) shows, there is relatively little change in the absolute value of the Debye length for the 
range of colloid diameters investigated. We therefore fix = 40 nm and vary the colloid 
size, and apply the same methodology as that outlined in sections HT1 and HVj in calculating 
the response for a = 400 — 2000 nm colloids to an external electric field. While an accurate 
determination of the melting transition is complicated by the small system size, melting 
is approximately determined by the splitting of the second peak in the radial distribution 
function g(r) u . We decrease the electric field to identify the weakest field at which the 
second peak in g{r) exhibits clear splitting [Fig. [9](a), inset], thus yielding d(E min ), the 
lattice spacing just prior to melting. The melting value of the lattice parameter d m is taken 
as 

d m = d(E mm ) + d{Emm 2 + 6E) (12) 

where 5E is the amount by which the field strength is varied between simulations, typically 
0.5 — 1.0 kVm~ l . The error in d m is then taken as d m — d(E min ). 

In terms of the particle diameter, d m grows considerably at small sizes, as the relatively 
larger Debye length leads to a minimum in the effective potential at relatively larger dis- 
tances. However, we recall that our analysis suggests a roughly a 3 dependence of e^p the 
strength of the dipolar attraction. This suggests that we might expect a stronger electric 
field to be required for relatively small colloids, and indeed we find that a higher field is 
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required to provide sufficient cohesive energy to hold the colloids together in the crystal-like 
'rafts', so at 100 kVm' 1 , a value we take as a reasonably accessible maximum field strength, 
the lattice spacing is still around d wo rs 1.6a for a = 400 nm, where d wo is the lattice spac- 
ing corresponding to a field strength of 100 kVm~ x . Conversely, at larger particle sizes, the 
Debye length is relatively small, so we find, even at melting, that the lattice constant only 
approaches around 1.2a. We note that, at high field strengths, the assumptions of Section 
UTI will ultimately break down. The results presented in Figs. [7] and [8] show that E < 100 
kVm~ l is reasonable, at least for a > 757 nm, supporting our conclusion of a maximum in 
tuneability around 1.2a. We show the values of the well depth of equation ( TTDT) at melting 
for the different experimental system in table HJ which indicates a trend towards deeper wells 
for larger particles, i.e. shorter ranges of the repulsive interaction relative to the particle 
size. 

Larger field strengths may be expected to yield greater tunability for the small particle 
sizes. Since the overall minimum lattice spacing is d — a, in principle smaller particles 
have greater tunability, however experimental observations reveal intense fluid turbulance 
disrupts crystal formation at higher field strengths. This is caused by large field gradients 
at the electrode edges that induce strong dielectrophoretic forces. In addition, higher field 
strengths can cause bubble formation by electrolysis and electrochemical degradation of the 
electrodes. 

The result is thus that intermediate particle sizes have the highest degree of tuneability, 
for a given maximum field strength. Here we have considered 100 kVm~ l , which leads to a 
maximum tuneability of d m /d\oo in the range 1000 nm < a <1500 nm. 

VII. CONCLUSIONS 

2D Colloidal crystalline 'rafts' with externally controllable properties have been modeled 
with Monte-Carlo simulations, and the resulting lattice constant, both experimental and sim- 
ulated, has been compared to the minimum in the effective interaction potential. In treating 
the effective colloidal interactions, we find that a straightforward, pairwise approach provides 
a reasonable description in describing the lattice parameter. The minimum in the effective 
interaction formed from combining the electrostatic repulsions and dipolar attractions is a 
useful means to approximate the lattice parameter, although close to melting fluctuations 
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lead to a larger value in the lattice parameter than this minimum. Conversely, at high field 
strengths, our treatment appears to be less accurate. Apparently, higher order terms which 
we have not accounted for may become important, such as non-linear Poisson-Boltzmann 
contributions to the electrostatic interactions, leading to a deviations from the Yukawa form 
[equation (pQ)}22 or limitations in our pairwise treatment of the dipolar attractions 23 . Other 
possibilities include inaccuracies in our assumptions for the effective colloid charge [equa- 
tion (j3J)], limitations of the 2D behaviour of the experimental system. Perhaps due to a 
cancellation of errors, at high field strengths, the minimum in the effective interaction can 
sometimes provide a more accurate value of the lattice parameter than the simulations. 

Another important assumption lies in our derivation on the electrostatic interactions. 
While their Yukawa-like form is well-accepted 2 ^, the values the effective colloid charge comes 
from equation Q. We note that the Debye length we determine seems rather constant across 
the four particle sizes considered, [Figj5] (a)], which suggests that out approach is at least 
consistent. However challenging it may be, a more quantitative measurement of the local 
ion concentration in the vicinity of the crystalline rafts is desirable and will be considered 
in the future. Furthermore, there may be some variation in the value of Z e ff. In fact this 
should affect all state points in a similar manner (for each particle size), and therefore we 
anticipate a similar outcome, but that a different value might be arrived at in the fitting of 
the Debye length, however, owing to the relatively short-ranged nature of the interaction, 
the effect of a different value for the effective charge is unlikely to impinge significantly the 
Debye length and should have a negligible effect on our main results. 

We have also assumed that a given system is described by one singleDebye length. In 
fact, in the counterion-dominated regime, the Debye length is in fact a function of colloid 
concentration 3 ^, an effect we have neglected. This might lead to a tightening at higher 
colloid concentration (i.e. high field strength), which may be expected to lead to an increase 
in deviation with the experimental results (Figs. [7]and[E]). Nevertheless, varying the Debye 
length as a function of local colloid concentration would be worth considering in the future. 

We have extended our approach a range of colloid size, to optimise the lattice tuneability 
for display applications. Tunable 2D crystal rafts can behave as tunable diffraction gratings 
capable of filtering white light into visible colours(author?) 11 . Similar tunable diffraction 
has been proposed for display devices 36 . The range of colours obtainable for a given geometry 
is governed by the lattice tunability. Fig 9b demonstrates that lattice tunability in our sys- 
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tern is maximised for particles of approximately 1 to 1.5 microns. Only higher field strengths 
applied to smaller particles can increase the lattice tunability. As already indicated, higher 
field strengths are not practical in the current experimental system, however in future ex- 
periments metal coated colloids would exhibit much stronger dipolar interactions, enabling 
a reduction in field strength for a given attraction. Regardless, our methodology illustrates 
using the well-developed machinery of effective colloidal interactions, as a means to model 
potentially complex interactions, useful for engineering purposes. 
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diameter 


Effective Yukawa 


KG 


Ionic 


C potential melting well depth 


(nm) 


contact potential (5e 




strength (mMol) 


(mV) 


(knT) 

\' V JD-*- j 


757±86 


4.60±0.05xl0 3 


19.0±1.0 


0.06±0.01 


-34.0 ± 1.4 


2.7±0.1 


965±44 


5.78±0.02xl0 3 


28.4±2.0 


0.08±0.01 


-37.1 ± 1.4 


2.9±0.1 


1390±30 


8.27±0.05 x 10 3 


32.3±3.0 


0.05±0.005 


-39.9 ± 2.0 


3.1±0.1 


2070±40 


1.17±0.01 x 10 4 


38.8±3.0 


0.04±0.005 


-43.6 ± 2.4 


3.4±0.1 



Table I: Fitted parameters of the Debye screening length used in the MC simulations for systems 
with different colloid diameters. Errors in the diameter are standard deviations in the data obtained 
found from electron microscopy measurements, while those in (3e, k<j, and the ionic strength are 
estimated from the fitting of d from MC data.£ potential measurements were made in a 0.01 mMol 
KC1 solution 
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Fi gure 1: The interaction potentials considered in this system. The Yukawa repulsion iiy«A;(r) 
models the electrostatic repulsions [dotted line, equation (JT])], the dipolar attractions are plotted 
as Udi p (r) dashed line, [equation 0]. These are combined to yield the total potential utot{r), solid 
gray line (equation HOp .This figure corresponds to a particle diameter of a = 757 nm and a field 
strength E = 30 kVm -1 

Figure 2: A schematic of the experimental set-up, side view (a) and plan view (b). The electric 
field is assumed to be constant across the observation region, as shown in (b). (b) depicts the 
quadrupolar electrode configuration (a,b),(c,d) 

Figure 3: Images of crystal rafts for the a = 945 nm system for applied fields of 29.0 and 80.0 
kVm~ l tox a and b respectively. Scale bars=10 [im.c snapshot from MC simulation for same system 
for a field of 30.0 jfeFm^and 2V~=144. 
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Figure 4: Change in the lattice constant for simulated systems of different sizes. Here a = 945 nm 
and a field strength E = 20 kVm~ l . Each run was repeated four times except ./V = 576. Error 
bars denote standard deviations over the results from different runs. Dashed line is solely a guide 
to the eye. 

Figure 5: (a) The Debye length k^ 1 fitted from Monte-Carlo simulation as a function of particle 
size. Dotted line denotes k _1 =40 nm taken for the simulations in section IVT1 Error bars denote 
the uncertainty in the Debye length resulting from the fitting in Figs. [7] and El (b) The contact 
potential of the Yukawa interaction (3e yu k [equation ([2])] as a function of particle size. We include 
a = 400 nm for which we take the Debye length k _1 =40 nm. 

The straight line is fit of the function f}£ yu k(cr) = Aa b where the fitted parameters are 

A = 11.1 ± 1.1 and b = 0.91 ± 0.01. 

Figure 7: Experimental results compared to MC simulation and the minimum in the effective 
potential, a 757 nm and b 945 nm diameter colloids. Filled squares are N = 72 while crossed 
are N = 144 from MC simulation. Circles are experimental data. Grey lines are the minimum of 
equation 1101 Error bars in the experimental data are standard deviations. 

Figure 8: Experimental results compared to MC simulation and the minimum in the effective 
potential, a 1390 nm and b 2070 nm diameter colloids. Filled squares are N = 72 while crossed 
are N = 144 from MC simulation. Circles are experimental data. Grey lines are the minimum of 
equation 1101 Error bars in the experimental data are standard deviations. 

Figure 9: (a) The lattice spacing around melting (dm, filled squares), and at at applied field of 
100 kVm~ l , as determined from MC simulation (crosses) and from Equation [TOj, as a function of 
particle size. Inset shows a radial distribution function close to melting (cr = 1000 nm, E = 11.0 
kVm^ 1 ). Arrow indicates split second peak. 

(b) the tunability ratio as a function of particle size. Connecting lines are guides to the eye. 

Figure 6: Radial distribution functions for an applied field of E=30.7 kVm~ 1 for the a = 2070 nm 
system. Circles are experimental data, line is MC simulation for A^=144. Arrows denote the first 
few peaks of a hexagonal lattice, d,V3d, 2d, V7d, 3d... 
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